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The dynamics of single bubble sonoluminescence (SBSL) 
strongly depends on the percentage of inert gas within the 
bubble. We propose a theory for this dependence, based 
on a combination of principles from sonochemistry and hy- 
drodynamic stability. The nitrogen and oxygen dissociation 
and subsequent reaction to water soluble gases implies that 
strongly forced air bubbles eventually consist of pure argon. 
Thus it is the partial argon (or any other inert gas) pressure 
which is relevant for stability. The theory provides quantita- 
tive explanations for many aspects of SBSL. 



Recent experiments |ij revealed that a single gas bub- 
ble levitated in a strong acoustic field P(t) = P a cos ujt 
can emit picosecond bursts of light, a phenomenon called 
single bubble sonoluminescence (SBSL). The phase and 
intensity of the light can be stable for hours. SBSL shows 
a sensitive dependence on the forcing pressure, the con- 
centration of the dissolved gas, and the liquid tempera- 
ture, among other parameters. 

Particularly puzzling is the dependence on the type of 
the gas dissolved in the liquid. Hiller et al. Q demon- 
strated that stable SBSL does not occur with pure nitro- 
gen, oxygen or nitrogen-oxygen mixtures. But a criti- 
cal concentration of argon gives stable SBSL. This paper 
presents an explanation for this dependence on the type 
of gas, based on combining principles from sonochemistry 
P,0 with hydrodynamic stability 

An important clue comes from experiments [HQ show- 
ing that the range of gas concentrations for stable SBSL 
in pure argon bubbles differs from that for air bubbles by 
two orders of magnitude. Stable SBSL in argon requires 
strong degassing of the liquid down to the tiny gas pres- 
sure Poo = 0.004Po7 which by Henry's law corresponds 
to a concentration Coo = 0.004co, where Pq = latm 
and Co is the saturation concentration. In contrast to 
argon, SBSL with air only requires degassing down to 
Pqo /P = 0.1 - 0.4 |,§. Lofstedt et al. @ estimated 
that diffusive equilibrium of the bubble with the sur- 
rounding dissolved gas requires Poo/Po ~ 10~ 3 , suggest- 
ing agreement with the experiments for argon, but not 
air. Indeed, detailed hydrodynamic stability calculations 
show quantitative agreement with argon data. The 
strong discrepancy led Lofstedt et al. to the conjecture 
that there is an "as yet unidentified mass ejection mech- 
anism" in air bubbles which "is the key to SL in a single 



bubble" . 

We suggest that this mechanism is chemical. The 
importance of chemical reactions has long been recog- 
nized in multibubble sonoluminescence (MBSL) in tran- 
sient cavitation clouds , since the pioneering work of 
Schultes and Gohr || found that aqueous solutions of ni- 
trogen produced nitric and nitrous acids when subjected 
to ultrasound. High temperatures generated by the bub- 
ble collapse are beyond the dissociation temperature of 
oxygen and nitrogen (« 9000K), leading to the forma- 
tion of O and N radicals which react with the H and O 
radicals formed from the dissociation of water vapor. Re- 
arrangement of the radicals will lead to the formation of 
NO, OH, NH, which eventually dissolve in water to form 
HNO2 and HNO3, among other products. 

Based on fits of SBSL spectra and hydrody- 

namic calculations p2|jr^ ], it is well accepted that inter- 
nal bubble temperatures in SBSL are even higher than 
in MBSL. Therefore, the same reactions as in MBSL will 
occur. The reaction products (NO2, NO, . . .) are pressed 
into the surrounding liquid, and are not recollected dur- 
ing the next bubble cycle, since their solubility in water is 
enormous. This chemical process deprives the gas in the 
bubble of its reactive components. Small amounts of N2 
and O2 that diffuse into the bubble during the expansion 
react and their dissociation products are expelled back 
into the surrounding liquids at the bubble collapse. The 
only gases that can remain within a SBSL bubble over 
many bubble cycles are those which at high tempera- 
tures do not react with the liquid vapor, i.e., inert gases. 
Hence, when air is dissolved in water, a strongly forced 
bubble is almost completely filled with argon. This argon 
rectification happens in SBSL but not in MBSL because 
it requires bubble stability over many oscillation cycles. 

In the following, we first present qualitative conse- 
quences of argon rectification, demonstrating that even 
at a crude level it resolves central problems of bubble 
stability. Then we proceed to make the argument more 
precise through quantitative calculation. 

If the bubble is filled essentially with argon, the hy- 
drodynamic stability of the bubble is determined by the 
partial pressure of argon p^' = £;Poo, where (— 0.01 for 
air) is the argon ratio of gas dissolved in the liquid. This 
fact immediately resolves the hundredfold difference be- 
tween the amount of degassing necessary for air versus ar- 
gon: As mentioned above, hydrodynamic stability calcu- 
lations Q demonstrate that stable sonoluminescence for 
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argon bubbles exists between pi OW er < P^/Po < Pupper, 
with [piower, Pupper] depending on the forcing pressure 
P a (figure @). E.g. at P a = 1.3atm, [p 

lower 3 Pupper \ ~ 

[0.002,0.004]. Since the stability window for air bubbles 
is set by the partial pressure of argon p££ — O.Olpoo, the 
total pressure of the air mixture must be in the range 
[lOOpiower, WOpupper] for stable SL. At P a — 1.3atm, this 
corresponds to 0.2 < Poo/Po < 0.4, in good agreement 
with experiments. One reason that air with its 1% argon 
is a particularly friendly gas for SL experiments is that 
this amount of degassing is easily achieved. 

FIG. 1. Phase diagram for pure argon bubbles in the 
Poo /Po versus P a /Po parameter space, from [^|, but now with 
experimental data included. Stable SL is only possible in a 
very small window of argon concentration. The experimental 
data points refer to observed stable SL (filled symbols) or sta- 
ble non-SL bubbles (open symbols) and are extracted (using 
the present theory) from refs. |7J (diamonds) and |H| (circles) 
and show good agreement with the theory. 

At this simple level, the theory makes several other ex- 
plicit predictions: When varying the percentage of argon 
in N2-Ar mixtures the range of Poo/Po where SL is sta- 
ble should vary like piower/Zl < Poo/Po < Pupper /fit, with 
the upper threshold, the lower threshold, and the range 
increasing with decreasing argon fraction Another 
consequence is that there is an argon ratio (f; « 0.003 
at P a — 1.3atm) for which stable SL should be possible 
without degassing. We caution that to achieve this it is 
necessary to rid the liquid of impurities to avoid sponta- 
neous cavitation. 

Also the different character of the transition to SBSL 
observed in air bubbles and in argon bubbles can be ex- 
plained by our theory: Hiller et al. || showed that for 
pure argon bubbles, the bubble mass increases smoothly 
and monotonically upon increasing the forcing pressure. 
In contrast Barber et al. || found that the transition to 
SBSL in air bubbles causes an abrupt decrease in the 
bubble mass. The difference between these two exper- 
iments is the presence or absence of sonochemical reac- 
tions: In air bubbles an abrupt behavior occurs when the 
forcing pressure corresponds to the onset of the dissocia- 
tion reaction. Below this threshold, the bubble contains 
air, and the ambient radius is determined by the diffusive 
stability of the mixture. Above the transition, when the 
molecular gases dissociate, the equilibrium radius is set 
by diffusive stability of pure argon bubbles resulting in 
a much different ambient radius [^). The transition to- 
wards SL is smooth for pure argon because a dissociation 
mechanism is absent. 

We now turn to a quantitative calculation of hydro- 
dynamic stability for gas mixtures with chemical reac- 
tions. The bubble radius R(t) is well described by the 
Rayleigh-Plesset equation 112], which we take to have 
the same form and parameters as in our earlier work || . 
The internal bubble pressure p(t) is assumed to obey a 
van der Waals equation of state. Now consider a bubble 



in water containing a mixture of a reactive gas (taken 
to be N2) and an inert gas, Ar. The total number of 
moles of gas in the bubble is N to t = 47ri?QP /(3G9o) = 
Nn 2 + N Ar , where Oo = 273K is the normal tempera- 
ture and G — 8.3143 J / (molK) is the gas constant. The 
argon fraction in the bubble is fb = N Ar /N to t, that of 
nitrogen 1 - f b = N N jN tot . If c Ar (r,t) and c N2 (r,t) are 
the mass-per- volume concentration fields of Ar and N2 in 
the liquid, respectively, the rate of change of the number 
of molecules of N2 and Ar in the bubble is given by 



N Ar = 
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Nn2 = ^ R2D »* drC N* \ r=R _ AN N2 exp . (2) 

Here, D Ar , Dn 2 , HAt and /iat 2 are the respective dif- 
fusion constants and molecular masses. The concentra- 
tion fields obey an advection diffusion equation, whose 
boundary conditions are set by the external concentra- 
tions c a (00, t) — = p^Cq I Pq (Henry's law) and by the 
partial gas pressures p a {t) in the bubble c a (R(t),t) — 
p a {R{t))c%/P , a = Ar, N 2 . The second term in @ 
represents the bubble's nitrogen loss by chemical reac- 
tions. The reaction rate depends on the temperature 
0(i) in the bubble. For simplicity, we assume that the 
reactions follow an Arrhenius law, with empirical pa- 
rameters appropriate for nitrogen dissociation (ref. Q ) : 
A = 6 ■ 10 19 (e /e) 2 - 5 (po/nN2)(Ro/R) 3 cm 3 /(mol S ) giv- 
ing the timescale of the reaction; 8* = 113000^ is the 
activation temperature and po the equilibrium gas den- 
sity. This reaction law is rather crude, as it neglects 
backward reactions as well as the kinetics of the expul- 
sion of reaction products; however, it is sufficient for this 
demonstrative calculation. 

The diffusive mass flux into the bubble can be calcu- 
lated explicitly using the fact that the diffusive timescale 
is much slower than the bubble oscillation period T 
|l7| , ^5| j7ptl . This reduces the diffusional problem to the 
calculation of weighted averages of the form (f) i — 

Jq f(t)R l (t)dt/ J Q T R l {t)dt, with the mass flux propor- 
tional to poo — (p) 4 for a pure gas ]lq| . Applying the 
same approximation to equations (|l|) and (|2|) gives 
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R m = max t i?(<). To close the equations, we need a 
model for the temperature dependence Q(t). The ac- 
tual temperature dependence is determined by compli- 
cated nonlinear processes operating during the collapse. 
As a simple model, we take the temperature to be uni- 
form within the bubble, and use the polytropic law 
0(t) = 6 {{R 3 - h 3 )/(R 3 {t) - h 3 ))'^ 1 with h the van 
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der Waals hard core radius and 7 the polytropic ex- 
ponent. The value of 7 depends on the rate of heat 
transport from the bubble, which is characterized by the 
Peclet number Pe = RRq/(Rk), where k is the ther- 
mal diffusivity of the gas. During the bubble expansion, 
Pe <C 1 (isothermal behavior, 7 = 1); during the bub- 
ble collapse, Pe ^> 1 (adiabatic compression, 7 = 5/3). 
Since both of these regimes occur during a single bubble 
cycle, it is necessary to use a model which interpolates 
between them. For definiteness, we follow Prosperetti 
p6| and use his calculated 7(Pe). We emphasize that al- 
though this treatment of the chemical reactions is crude, 
the central results discussed below are robust. 

FIG. 2. Phase diagram for air at Poo/Po = 0.2 in the Ro 
- P a space. The arrows denote whether the ambient radius 
grows or shrinks at this parameter value. Curve A denotes 
the equilibrium for an air bubble, on curve C the bubble only 
contains argon. The intermediate curve B necessarily exists 
because of the topology of the diagram, and represents an 
additional stable equilibrium. Above and right of the thin 
line, the gas temperature exceeds the nitrogen dissociation 
threshold of about 9000K. 

With these approximations, the equilibrium states 
(satisfying AN Ar = ANjy 2 = 0) can be computed as 
a function of (Na t , Nn 2 ), or equivalcntly as a function 
of (£{,,i?o). Fig. H shows the equilibrium radii Rq in 
the Rq — P a plane for air at Poo/Po = 0.2. For small 
forcing the temperatures are not high enough to initiate 
chemical reactions, so that the equilibrium curve corre- 
sponds to the classical prediction by Eller and Flynn |T^] 
for this gas concentration. This equilibrium is unsta- 
ble: The bubble either shrinks or grows by rectified dif- 
fusion; experiments show that a growing bubble 
eventually runs into a shape instability where microbub- 
bles pinch off and make the bubble dance because of the 
recoil In the opposite limit of high forcing (curve 
C), the reactions burn off all the N2, so that the bub- 
ble contains pure argon; this equilibrium corresponds to 
the (stable) Eller-Flynn equilibrium at the argon partial 
pressure p£/P = OMp^/Po = 0.002. 

Figure ^| displays a regime of shrinking bubbles at high 
forcing pressures (left of curve C) and an adjacent region 
of growing bubbles (right of curve A) . This necessitates 
the existence of an additional equilibrium at intermediate 
forcing pressures, curve B in figure |^, for which growth 
by rectified diffusion and loss by reactions balance. This 
additional equilibrium occurs close to the point of nitro- 
gen dissociation, and turns out to be stable; the argon 
fraction for this equilibrium is slightly larger than the 
fraction in the liquid (for not too strong forcing). - 
The only feature of Fig. 2 that depends on the details of 
temperature and chemical reactions is the exact position 
of the nitrogen dissociation threshold and thus the exact 
position of curve B. 



FIG. 3. The equilibrium fraction of argon in the bubble 
as a function of P a and Ro. 

Figure |^ plots the equilibrium composition Q , given 
by AN Ar /AN N2 = - Weakly forced bubbles 

have Q * £/, thus Poo/Po — 0.20 is relevant for stabil- 
ity. Strongly forced bubbles have w 1 £/, thus 
Pod/Po — 0.002 is the relevant quantity. The transition 
between these regimes is abrupt, and occurs when the 
bubble temperature surpasses the dissociation tempera- 
ture (» 9000^ for N 2 ). 

What happens for even lower argon concentration £j < 
0.01 in the dissolved gas? For these low concentrations 
the equilibria curves A and B hardly depend on . This 
holds even in the limit — > of pure nitrogen bubbles. 
Therefore, our theory predicts that there is a parameter 
regime of forcing pressures where stable N2 bubbles exist. 
The equilibrium curve C, of course, does depend on and 
for decreasing argon concentration p^ = £;Poo it moves 
further and further to the right, allowing diffusively sta- 
ble SL bubbles only for larger and larger forcing. Finally, 
spherical instability will destroy these bubbles ||. 

Having calculated the above phase diagrams, we come 
back to a comparison to experiments. For an air bub- 
ble, upon increasing the forcing pressure, the bubble 
will eventually encounter the stable equilibrium curve B; 
tracking curve B on further increasing the forcing pres- 
sure leads to an abrupt decrease in the bubble size. The 
experimentally observed abrupt decrease in the ambient 
radius reflects the sharpness of the slope of curve B. Even- 
tually, at high enough forcing pressure, the bubble will 
track curve C, the ambient radius now increasing with 
forcing pressure. 

The most extensive experimental evidence in support 
of these calculations is Holt and Gaitan's recent de- 
tailed measurement of phase diagrams in the Rq — P a 
plane. At low forcing pressure, they observe bubble equi- 
libria in agreement with classical calculations using the 
pressure head (poo/Po=0.2) applied in the experiments. 
At very high forcing pressure, the bubble equilibria agree 
with classical calculations, but only when p^ / P = 0.002 
is used in the equations. Between these two regimes, Holt 
and Gaitan find the additional stable equilibria (curve B) 
of non-sonoluminescing bubbles and an adjacent region 
(at P a ~ 1.2 — 1.3atm) where bubbles dissolve, very simi- 
lar to that between curves B and C in Fig. ||. In addition, 
Holt and Gaitan also find the the size of the dissolution 
region decreases with increasing Poo/Pq, which also fol- 
lows directly from the theory. 

To summarize, the combination of the principles of 
sonochemistry and hydrodynamic stability leads to a con- 
sistent picture of the stability of sonoluminescing bubbles 
for gas mixtures and makes many predictions. Our cen- 
tral statement is that it is only p^/Pq which is rele- 
vant for SL stability in the high P a regime. We included 
all available experimental data (with sufficient informa- 
tion on P a and p^) in our theoretical phase diagram 
H figure 1 to demonstrate the good agreement. Fi- 
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nally, we suggest to measure the concentrations of the 
reaction products as a function of time for SBSL, as al- 
ready done in MBSL Q . Nitrous acid production would 
lead to a decrease in pH. For an estimate of an upper 
bound to the production rate we assume that reactions 
at the collapse burn off all the nitrogen that diffuses into 
the bubble during its expansion. This amount of gas 
is estimated in ref. as AA^v 2 = 2irD]y 2 c^R m T/ /j,n 2 
per cycle. With typical values of R m — IORq for the 
maximal radius, Rq = 5/im, D^ 2 = 2 ■ 10 _9 m 2 /s, 
« 0.20c^ 2 , c^ 2 = 0.02fc 5 /m 3 , and T = 37/is one 
obtains A7Vjv 2 ~ 3 • 10~ 18 mol per cycle or ~ 3 • 1CU 10 
mol of N2 per hour converted to reaction products. This 
results in a small but detectable pH decrease. 
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